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ABSTRACT 

Context. In Unification Models, Active Galactic Nuclei (AGN) are believed to be surrounded by an axisymmetric structure of dust and gas, 

which greatly influences their observed properties according to the direction from which they are observed. 

Aims. The main aim of this work is to constrain the properties of this obscuring material using X-Ray observations. 

Methods. The distribution of column densities observed by Chandra in the Chandra Deep Field South is used to determine geometrical 
constraints for already proposed torus models. 

Results. It is found that the best torus model is given by a classical 'donut shape' with an exponential angular dependency of the density 
profile. The opening angle is strongly constrained by the observed column densities. Other proposed torus models are clearly rejected by the 
X-Ray observations. 
Conclusions. 
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1. Introduction 

A unified scheme to explain a variety of Active Galactic Nuclei 
(AGNs) has been proposed by Antonucci ( 1993) to explain the 
differences observed in Seyfert galaxies, a class of local ac- 
tive galaxies. According to this Unification Model, the active 
nucleus is surrounded by a toroidal structure composed of gas 
and dust, which determines dramatic spectral differences de- 
pending on our line of sight towards the central source. In the 
X-Ray domain, if an AGN has a neutral Hydrogen column den- 
sity {Nh) in the line of sight smaller (when viewed face-on) or 
larger (when viewed edge-on) than 10 22 cirT 2 , then the object 
is classified as Type I or Type II, respectively. 

The dusty toroidal structure that surrounds the central 
source absorbs the optical, ultraviolet and soft X-Ray emission 
from the active nucleus. However, hard X-Ray photons (> 2 
keV) can partially escape from this cold material. Therefore, 
from X-Ray observations it is possible to measure the column 
density in the line of sight, and put tight constraints on the prop- 
erties of the obscuring medium. 

It is commonly believed that this X-Ray emission is pro- 
duced in a 'corona' near the super-massive black hole, where 
low energy photons from the accretion disk are reprocessed 
by energetic electrons (either mildly relativistic thermal elec- 
trons or highly relativistic non-thermal electrons) via Inverse 
Compton Scattering. The main result of this process is an ob- 



served spectral energy distribution described by a power law, 
with a typical slope F « 1.8 (Turner et al. 119971 Tozzi et al. 
2006), and an exponential cutoff energy at w 300 keV (Matt et 
al.[1999). The radiation can also be reflected and/or scattered 
depending on the circumnuclear material distribution causing 
the overall spectral shape to become flatter, F w 1 .7, with an 
apparent smaller energy cutoff (see Svensson ll996l) . 

As hard X-Ray energy photons can penetrate the surround- 
ing material, and eventually escape from the AGN host, they 
can be detected with the present generation of X-Ray tele- 
scopes. But if the Hydrogen column density along the line of 
sight becomes larger than the inverse of the photoelectric Cross 



Section o-^y 



10 24 cirT 2 (see Matt 2002 ), then the medium 



becomes Compton Thick (t = Nhct ~ 1). In this case, the only 
observable emission component is coming from photons scat- 
tered and/or reflected by the circumnuclear media (Wilman & 
Fabian [r999l ). 

Since Compton Thin AGNs (those with N H < 10 24 cirT 2 ) 
are ubiquitous emitters in the 2 - 10 keV hard X-Ray band, deep 
X-Ray surveys have proven to be the best way to estimate the 
number density and evolution of active galaxies in the Universe 
(Mushotzky 2004). Unfortunately, these X-Ray surveys are al- 
ways flux limited (nowadays ^"uikev ~ 10~ 16 er E sec 1 cm -2 ), 
and therefore biased against faint sources. In fact, the popu- 
lation of high redshift, low luminosity, and highly obscured 
AGNs is still not well determined. This problem translates into 
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a poorly determined faint end of the luminosity function, im- 
plying uncertainties in the contribution of obscured AGNs to 
the Cosmic X-Ray Background (CXRB), which is known to 
require of a sizable population of obscured sources in order to 
explain its hard spectral energy distribution. 

The aim of the present work is to constrain different 
torus geometries that can reproduce the observed distribution 
of Hydrogen column densities (Nh) in AGNs observed by 
Chandra in the Chandra Deep Field South (CDF-S). 

Section 2 describes the intrinsic Nh distribution found by 
Tozzi et al. (f2006) which was used in our analysis. In Section 
3 we explain the theoretical modelling, while in Section 4 
we present our models and results. In Section 5 our results 
are discussed and we present the conclusions in Section 6. 
Throughout this work we use: Q„, = 0.3, £2a = 0.7 and H = 70 
km sec^ 1 Mpc" 1 (Spergel et al. 12005V 

2. Deep X-Ray Observations 

2.1. AGNs in the Chandra Deep Field-South 

In the past few years, several deep Chandra images of the ex- 
tragalactic sky have been obtained, with the 2 Ms Chandra 
Deep Field-North (CDF-N; Alexander et al. l2003l and 1 Ms 
Chandra Deep Field-South (CDF-S; Giacconi et al. l2002l be- 
ing the two deepest (F^J keV » 1.4 • 10~ 16 erg s cm 2 for CDF- 
N and a factor 2 shallower for the CDF-S). Optical identifi- 
cations have been obtained for most X-Ray sources in these 
fields. Barger et al. (2003 ) presented multicolour imaging data 
for all 503 X-Ray point sources in the CDF-N and spectro- 
scopic redshifts for ~ 56% of the sample. Zheng et al (2004) 
presented accurate photometric redshifts for 342 sources in the 
CDF-S using the large multicolour data set available for this 
field. This corresponds to ~ 99% of the sample and includes 
173 sources with reliable spectroscopic redshifts from Szokoly 
et. al (120041 

Given the photo-z completeness of ~ 99% for the the X- 
Ray sources in the CDF-S, this field becomes the perfect choice 
for our study of the torus properties. Redshifts are necessary in 
order to compute the observed Nh (see below) as well as the 
source X-Ray luminosity. 

Following the sample definition for the CDF-S found in 
Tozzi et al. ( 120061 ), out of the 347 X-Ray detected sources (346 
sources presented by Giacconi et al. ( 2002 ) and 1 extra source 
presented by Szokoly et al. (2004)), 7 sources are stars and 4 
have no redshift information. Further 15 sources have X-Ray 
luminosities < 10 41 erg s _1 cirT 2 , a luminosity range where 
emission from star forming galaxies becomes significant, and 
are therefore discarded. Finally, 14 Compton Thick AGN were 
also eliminated, as discussed in Section 4. The final sample 
contains 307 X-Ray point sources at z > that will be used for 
our study. 

2.2. Hydrogen Column Densities 

The distribution of Nh values for the sources in the CDF-S 
were derived by Tozzi et. al (2006). They performed direct 
X-Ray spectral fitting to the observed Chandra data using the 



redshift information published by Zheng et. al (2004). Using 
a subsample with the 82 brightest sources, Tozzi et. al ([2006) 
determined a weighted mean spectral index of T ~ 1 .75 which 
resulted independent of the values of the derived Nh- By fix- 
ing F = 1.8 the redshifted Nh was derived for the rest of the 
sample. The distribution was then corrected for incompleteness 
in the luminosity and redshift parameter space using the sur- 
vey effective area and the luminosity function of Ueda et al. 
(2003). Ueda's luminosity function was determined using 247 
AGN detected in the 2-10 keV X-Ray band from a compi- 
lation of ASCA, HEAO - 1 and Chandra observations. The 
luminosity function shows a luminosity dependent density evo- 
lution, in which the low luminosity AGN population peaks at a 
lower redshift than the high luminosity sources. This is consis- 
tent with optical quasar observations in which this population 
peaks at redshift ~ 2 (Boyle et al. l20001 l. 

The final corrected Nh distribution, seen in Figure 16 of 
Tozzi et. al (2006), will be used in the following sections to 
constrain the torus parameters. 

3. Model General Properties 

We first assume that the X-Ray emission in AGNs comes from 
the nearest region to the super-massive black hole. This re- 
gion has a typical spatial scale of Rq ~ 10 - 100 Rs (where 
R$ = GMbh/c 2 is the Schwarzschild Radius), implying phys- 
ical sizes of Rc ~ 10~ 5 - 10~ 3 pc for a typical super-massive 
black hole with M BH ~ 10 7 - 10 8 M Q . Recent observations of 
NGC 1068 using mid-infrared interferometry, suggest a torus 
about 3.4 pc diameter (Jaffe et al. 120041 ). while near-IR rever- 
beration mapping of nearby Seyfert galaxies estimates an inner 
limit at ~ 10~ 2 - 10~' pc (Suganuma et al. 12006V Therefore, 
for practical purposes we can approximate the X-Ray source 
as a point-like emitting region, implying an easy treatment for 
the X-Ray radiation that will provide the total column density 
along the line of sight. 

3.1. Torus Properties 

We consider four different torus models (presented in Section 
H]i to obtain synthetic distributions of column densities which 
are then compared with the observed distribution from Tozzi 
et. al (I2006V 

We model the obscuring region as a simple axisymmetri- 
cal matter distribution surrounding the accretion disk. Given 
the symmetry of the problem, random lines of sight with polar 
angles between and n/2 (weighted by the angle differential 
area), are used to calculate the optical depth for each random 
direction, where the the polar angle <p is defined as the angle 
subtended between the line of sight and the torus equatorial 
plane. 

We assume a solar abundance for the torus, a photoelectric 
cross section for absorption given by Morrison & McCammon 
(11983V a standard Galactic Gas-to-Dust ratio N H /E(B - V) = 
5.8 • 10 21 cm -2 given by Bohlin et al. ( 1978) and a solar neigh- 
bourhood value Ry = 3.1 (Schultz & Weimer I"l975l) . These 
assumptions relate the Hydrogen column density Nh with the 
optical extinction Ay, and the optical depth Ty (Ay - 1 .09tv), 
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NhW) = <r _1 T = J p(r,4>)dr 

= 5.8 • 10 21 E(B - V) cm" 2 = 1.9 ■ 10 21 A V cm" 2 (1) 
= 2.0 • 10 21 t> cm" 2 = 2.0 • 10 21 T Vi N T cm" 2 

where, if the medium is continuous, p(r, 0) is the matter den- 
sity distribution and <x the effective absorption cross section. 
Otherwise, if the medium is clumpy, T Vi represents the optical 
depth for one cloud and Nj the average number of clouds along 
the line of sight. It is easy to show that if the density distribu- 
tion p has no angular dependency (i.e., p = p{r)), then the radial 
dependency will be hidden by the assumed maximum column 
density NHmax along the line of sight (see SectionlU). 

Other model parameters, such as the inner radius, exter- 
nal radius, density in the equatorial plane, and/or cloud opti- 
cal depth, are also constrained by NHmax- In this work we fix 
the maximum value for the column density at Nnmax = 10 25 
cm" 2 , according to X-Ray observations of Type II AGNs in 
the 0.1 - 100 keV band using the BeppoSAX X-Ray satellite 
(Maiolino et al. [1998b . 

Due to the degeneracy of the model parameters, we use the 
minimum number possible of free parameters to describe the 
torus distribution. The chosen parameterisations do not give 
physical sizes of the distributions explicitly, but instead re- 
duced parameters such as Rim I Rom (where R out is the exter- 
nal radius and /?,„, is the inner radius for the torus) are used. 
Ri„t can be estimated by the 'evaporation radius' for graphites 
(at T as 1500 K) defined by assuming thermal equilibrium be- 
tween the ultraviolet incident radiation and the rate for reemis- 
sion from the dust. Assuming a gr = 0.05 pm as a typical grain 
size, and following Barvainis ( 1987]), 

R int « 0.4L^ /2 45 pc, (2) 

where Luv.as is the ultraviolet (UV) luminosity per unit of 10 45 
erg sec" 1 . Assuming this value as a lower limit for the distance, 
we can estimate the external radius R out , and other model pa- 
rameters using our torus model results (see estimations in Table 

ID- 

4. Modelling the Torus 

The large degeneracy of model parameters and the uncertain- 
ties given by the lack of faint sources in the Nh distribution 
determined by Tozzi et al. ([2006), require us to keep the num- 
ber of free parameters within our torus models to a minimum. 
As such, we do not allow the parameters to change with lumi- 
nosity or redshift. 

Barger et al. ([2005 ) presented a complete sample of local 
X-Ray selected AGN (z < 1 .2), in which a luminosity depen- 
dence for the fraction of Type II over Type I AGNs is seen, as 
Ueda et al. ([2003 ) proposed earlier. This implies a possible de- 
pendency of the covering factor of the dusty torus on luminos- 
ity, and therefore a change in the geometrical torus distribution. 
Hence, our results can be understood as an average geometrical 
distribution for the entire AGN population, given by a simple 



axially symmetric geometry, with no dependency on luminos- 
ity or redshift. A discussion on the luminosity dependency can 
be found in Section[5] 

To determine the torus parameters which best fit the ob- 
served column densities, we use the x 1 test assuming poisso- 
nian error bars for the observed data. Our own experience deter- 
mining the intrinsic Nh distribution based on Hardness Ratio 
measurements suggests that this method introduces large un- 
certainties in the inferred Nh values, which dominate over the 
poissonian errors from a population of sources. On the other 
hand, the spectral treatment introduced by Tozzi et al. ((2006) 
gives more reliable values for the intrinsic properties of the ob- 
scuring region present in AGNs. 

Tozzi et al. d2006) found 14 sources showing a reflection- 
dominated spectrum. These Compton Thick sources were not 
used in our analysis because of the large uncertainties in the 
number counts of this population. In fact, we have not used 
the last bin in the intrinsic column density distribution shown 
in Figure 16 of Tozzi et al. 2006, which includes objects with 
Nh > 10 24 cm" 2 and corresponds to a lower limit only. 

The models presented in the remaining of this section were 
motivated by previous work which were carried out to explain 
a variety of observational properties of AGNs (Treister et al. 
120041 Pier & Krolik [19921 [TW3l Granato & Danese [19941 
Elitzur et al. 2003). A summary of the results from the mod- 
elling is presented in Table Q] 

4.1. Model 1 

Following the torus distribution proposed by Treister et al. 
( 2004), we modelled the geometry given by Figure [T^ to find 
different optical depths as a function of the line of sight, (f>. In 
this case the matter density is not distributed as in an homoge- 
neous medium, but has the following dependence with <p and 
r. 

p(n0)=Pe q (^f^ |SmWI 0) 

where p eq is the inner density (at r mt = R m -R t ) in the equatorial 
plane (constrained by NHmax, see Equation In this model, 
the free parameters are R m /R, (see Figure QJ), y, and p. The 
analytical equations which describe the column density as a 
function of <f> are given by: 

If < 4> < sin" 1 (R,/R m ) and /? * 1 then, 

N H = N Hmax e-v^W x 

, xi-fl / , \H>' 

cos(0)+ y/(R,/R,„) 2 -sm(4>) 2 \ -lcos(<W- \(R,/R m ) 2 -mn((l>) 2 \ 

( 1 +R, /R„, ) 1 ' -( 1 -R, IR,„ ) 1 ^ 

If < (f) < sin" 1 (R,/R m ) and = 1 then, 

Nh = N Hma x e" y|si " wl X 

ln((cos(0)+ V(«,/«™) 2 -sin(<» 2 )/(cos(<»- yJ(R,/R m ) 2 -sm(^ 
hm+R,/R,„)/(l-R,/R,„)) 

If sm- l (R t /R m )<4><n/2 then, N H = N Hmm 

(4) 

The maximum and minimum column density along the line 
of sight is through the equatorial plane and at the poles, respec- 
tively. 
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If P + 1 then 

Num. 



,1- 



a R m 

P e( i r iat ] _ p 



1 + 



R, 



If P = 1 then, AW« = Peq In 



^n. / \ R, 

1 + fl,//? m 



t-jS 



(5) 



1 -R t /R„ 



This model gives a good fit to the data, with v . = 1.30. 
Results adopting /? = 0.0 are presented in Figure [2] and show a 
well constrained region of the R m /R t - y parameter space. The 
best fit values are spread between 1.07 < R m /Ri < 1.15 and 
7.0 < y < 9.0 at 3 <x confidence levels. The same results are 
obtained when a value p = 2.0 is adopted. Adopting f3 = 1.0 
and 3.0 results in a small decrease in the value of y, but does 
not introduce further changes. Note however, that p modifies 
the value of Njimax through Equation[5] 

From the family of parameters found, we show in figure 
[3] some examples of synthetic Nh distributions and compared 
them with the one determined by Tozzi et al. (120061 . It is seen 
that the parameter R,„/Ri mostly affects the distribution of Type 
I sources (i.e., for Nh < 10 22 cm -2 ), while y changes the Type 
II/I fraction. 

A sketch of the torus geometrical distribution, choosing the 
best fit parameters, is given in Figure HI where the 'donut' 
structure of the model is clearly seen. Under these assumed 
parameters, this model predicts an intrinsic fraction of sources 
with column densities > 10 24 ctrT 2 of ~ 27%. 

4.2. Model 2 

This model is proposed following previous studies by Pier & 
Krolik ( 1992, 1993 ) based on the IR reemission coming from 
the torus. We use their geometrical distribution (Figure QJ)) to 
find the dependency of the column density with the line of 
sight. 

This model assumes a homogeneous density with hja and 
b/a as free parameters, where h, b, and a are as shown in Figure 
QJ). The dependency of the column density with the line of sight 
(0 < (p < 7i/2) is given by: 



If O<0<tan- 1 (A) th en, N H = N Hmax ^^ 



If tan -1 (A) < < tan -1 (A) then, 



Nh = Nh 



cos(tan~'(A/2fc)) / /, 



b/a-\ 



(Asin^)- 1 -cos(^)- 1 ) 



If tan" 1 (A) < (j) < tt/2 then, N H = N Hmin 

where the maximum column density is: 

b — a 

N Hmax — Peq 



(6) 



(7) 



cos tan 



Comparing the synthetic distributions of column densities 
with the observed one, we found the results shown in Figure 
[2] The parameters are not well constrained and imply 2.4 < 
hja < 3.3 and 250 < b/a < 600 at 3 <x confidence levels. The 



model also shows an extreme radial-to-vertical thickness ratio 
(b/h * 75 - 250). 

Synthetic Nh histograms are shown in Figure [3] It is seen 
that the parameter b/a systematically overestimates the number 
of sources at around Nh = 10 225 cm' 2 , while hja modifies the 
first bin of unobscured sources. 

Clearly, this model is limited by its constant density distri- 
bution. This translates into a large b/a, needed to reproduce the 
wide range of observed column densities (since Ng °c As), and 
a very sensitive distribution in the innermost region of the torus 
(see Figure|Ui. Also, the model does not predict a large fraction 
of highly obscured sources, implying that it cannot explain the 
CXRB or the number of Compton Thick sources observed in 
the CDF-S (~ 10%; Tozzi et al. [2006b . The typical intrinsic 
fraction of Compton Thick sources predicted by the model is 
only ~ 3% which are produced by a narrow angular section 
near the equatorial plane able to account for large column den- 
sities (see Figure[4|i. 

When we compare our analysis with the results by Pier & 
Krolik (1992), we find further disagreements. Typical dimen- 
sions deduced by them, as restricted by IR torus reemission, 
are of parsec-scale, meanwhile we obtain a torus with typical 
scales of ~ 200 pc assuming Ljjy = 10 45 erg sec -1 . 

We conclude that this torus model can be rejected due to 
the numerous caveats already discussed, such as typically large 
outer radius, which is not seen in local AGNs (Jaffe et al. 2004; 
Prieto et al 120041 ). the sensibility of the geometrical distribu- 
tion in the inner regions, and the small number of predicted 
Compton Thick sources. 

4.3. Model 3 

This model has the spatial density distribution shown in Figure 
QJ and is based on the work presented by Granato & Danese 
( 1994). The density profile has an exponential dependency with 
the polar angle parameterised by a factor y (when <p < <p c ), and 
a power law dependancy with the radius parameterised by the 
index /3: 



p(r,(f>) =p eq (^yj 



,-yMn(< 



(8) 



where p e q is the inner torus density at the equatorial plane. The 
free parameters of the model are the opening angle C , y and 
(3. The equation which describes the dependence of the column 
density with the line of sight is: 



(9) 



If <0<0 C then, NH=N Hmax e-^<M 
If <p c < </> < tt/2 then, Nh = N H min 

where Nnmin and Nnmax are fixed parameters at 10 20 and 10 25 
cm -2 , respectively. The minimum value fixes Nh at the poles 
and the maximum value fixes Nh at the equatorial plane. 



If p + 1 then, N H max = Peq >f nl 



/ 1-/8 1-/3 \ 
r — r 
out mt 



If p = 1 then, N H max = Peq r mt In 



1-p 



('"out 
rmt 



(10) 
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Model 


2 

A. vjnin 


Parameter 1 


Parameter 2 


Parameter 3 


Routlpc] 


Peq[cm~ 3 ] 


1 


1.30 


1.07 < RJR, < 1.15 


7.0 < y < 9.0 


0.0 < p < 3.0 * 


5.7 - 12 


2.8 - 6.0 ■ 10 s f 


2 


8.13 


2.4 < h/a < 3.3 


250 < b/a < 600 


X 


100 - 240 


1.3-3.2- 10 4 


3 


5.65 


1.0 < <p c rod < 1.2 


8.0 < y < 10.0 


free 


1.0* 


5.4 • 10 6 1 


4 


2.65 


0.38 <crrad< 0.44 


T Vi < 0.8 


X 


free 





Notes .- 

* indicate a fixed parameter. 

T these parameters can be modified taking into account radial dependencies for the model, i.e. p + 0. 

* obtained assuming y = 9.0. 

Table 1. Summary of the estimated parameters for the models described in Section|4] Ranges are obtained from 3 cr confidence 
limits (see contour plots in the Figure|2]l. The last two columns are deduced considering a UV luminosity of 10 45 erg sec -1 and 
Equation [2] No radial dependency for the density distribution (fj = 0.0) is assumed in the models. In Model 3 we use a fixed 
parameter R out — 1 pc to estimate a density at the equatorial plane, although this is a free parameter. Model 4 is not shown in the 
last columns due to the different cloud treatment explained in Subsection l4.4l 



Searching for best fit parameters using the x 1 test, we ob- 
tain the results shown in Figure [2] This model gives a poor 
fit to the data ixl min = 5.65), and is not able to reproduce 
the observed Nh distribution between the first bin and those at 
% 10 22 crrT 2 (Figure|3). The sharp gap seen at Nh ^ \0 2l cmT 2 
in the synthetic distributions corresponds to the transition pro- 
duced when the line of sight no longer hits the torus (0 > <p c ). 
This dichotomy also affects the x 2 distribution where multiple 
minima are seen. Parameters between 57° < <p c < 69° and 
8.0 < y < 10.0 are found. 

Simulations of the intrinsic Nh distributions are shown in 
Figure[3] Changes in the value of C affect the distribution over 
Type I sources only (i.e., for Nh < 10 22 cm -2 ), while y affects 
the number of sources at Nh ~ 10 21 cm~ 2 and changes the frac- 
tion of Type II/I. The estimated intrinsic fraction of Compton 
Thick sources, using the best fit parameters, is about 26%. 

When we compare our results with the work by Granato 
& Danese (1994), we find that our analysis can constrain <p c 
(except for large y values), while Granato & Danese (1994) 
found that the spectral energy distribution for the reemission 
does not depend strongly on this parameter. On the other hand, 
we cannot constrain the outer radius, since this scale parameter 
is contained in Nnmax (EquationfTOli, and it could be easily ad- 
justed to adopt any value (e.g. > 100 pc, as found by Granato 
& Danese [l994i ). 

4.4. Model 4 

Many previous theoretical models trying to reproduce the IR 
reemission coming from the obscuring torus were based on 
continuous density distributions (as those presented in Models 
1 to 3). Elitzur et al. ( 120031 and Nenkova et al. ([2002b have 
developed the formalism to handle dust dumpiness and their 
results indicate that the inclusion of thick clouds may resolve 
the difficulties encountered by previous theoretical efforts try- 
ing to reproduce the observed IR spectral energy distributions 
of local AGNs. They have shown that some problems are nat- 
urally resolved if a clumpy media with ~ 5 — 10 clouds along 
radial rays (each one with t v . > 40) is considered. The fun- 
damental difference between clumpy and continuous density 
distributions is that radiation can propagate freely between dif- 



ferent regions of an optically thick medium when it is clumpy, 
but not otherwise. 

Following Elitzur et al. (2003) we define the model given 
by Figure[TJl, considering a spatially random cloud distribution 
per unit length given by: 

n T (r,4>) = n eq ^fJ (11) 

where n eq is the inner number of clouds by unit length in the 
equatorial plane, cr the half opening angle for the angular dis- 
tribution, and ft the power law index describing the radial de- 
pendency. Integrating the last equation along the line of sight, 
we obtain the average number of clouds, Nt as a function of 
the line of sight, parameterised by the angle <p: 

N T (<f>) = N eq e-?^ (12) 

where N eq is the average number of clouds along the equatorial 
plane: 



If ft+\ then, N eq = n eq -^-(r 1 J-r^) 

(r (B) 
If /J = 1 then, N eq = n eq r int In — 

V r int I 

Assuming ry = NrTv n and considering EquationQ] we can 
estimate the Hydrogen column density dependency on the line 
of sight, as: 

N H (4>) = NH m a X e-* 2/cT2 (14) 

The fixed parameter Nnmax constrains the optical depth, the 
number of clouds along the line of sight at the equatorial plane, 
and the parameter /3: 

N Hmax = 2.0 ■ 10 21 jV eq T Vi (15) 

If a certain value T Vi is assumed for the clumps, the only 
remaining free parameter for this distribution is cr. However, 
restricting Nt to represent a discrete number of clouds and 
assuming a poissonian probability for the number of clouds 
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along the line of sight, it is also possible to constrain the op- 
tical depth per cloud, tv,. In the case when the optical depth 
tends to zero but Nj tend to infinity, the continuous case is re- 
covered, where column densities are allowed to have any value. 
Note that Nenkova et al. (2002 ) assumed Ty f = 40 to reproduce 
the IR reemission from the torus, implying a minimum col- 
umn density of 8.0 • 10 22 cm 2 . Clearly, smaller optical depths 
per cloud are necessary to explain the observed distribution of 
Nh- A work around this difficulty which would still be consis- 
tent with a large Ty,- per cloud, is to assume that absorption by 
overdensities of diffuse material located in the outer regions or 
within the torus might be responsible for the smaller observed 
column densities. 

The best-fit values for cr and Ty; are shown in Figure 
12 Multiple minima in the x 1 distribution are found due to 
the cloud number discretization and the selected Nh bin-size. 
Values between 22° < cr < 25° and ry; < 0.8 are preferred. 
Parameter cr changes the number of objects at Nh = 10 20 cm" 2 
and at Nh = 10 25 cm" 2 , without altering significantly the inter- 
mediate values where a increasing slope in the synthetic distri- 
butions is seen. On the other hand, the optical depth per cloud 
creates a dichotomy between Type II and Type I objects, where 
the minimum value of Nh is given by the optical depth of a 
single cloud. 

From Figure 2J it is clear that this model gives the largest 
number of Compton Thick sources produced by lines of sights 
near the equatorial plane (approximately a fraction of ~ 58% 
using best fit parameters). On the other hand, if a lower Nnmax 
of 10 24 cm" 2 is adopted, the gaussian angular dependency gen- 
erates large discrepancies with the observed column densities. 

When we compare our parameters with the results coming 
from the reemission treatment given by Elitzur et al. ([3003 ), 
we found that our typical Gaussian angular distributions are 
smaller than Elitzur's values. In fact, they find a valid range 
of cr = 45° ± 15°. Also, their estimations for the number of 
clouds along the equatorial plane (N eq ~ 5 - 10) correspond to 
a maximum column density of ~ 10 24 cm" 2 (see Equation [1), 
one order of magnitude smaller than our assumptions. In our 
models, since NHmax - 10 25 cm" 2 , the value of N eq depends on 
the assumed Tyi- 

Note that since Nh Nt, this clumpy model requires a 
wide range for the number of clumps (at least ~ 3 orders of 
magnitude) to describe the wide range of observed A^s. 

5. Discussion 

5.1. Torus Properties as a Function of Luminosity 

The work of Barger et al. (120051 . based on a complete AGN 
sample at z < 1 .2, shows a dependency of the fraction of Type 
II/I AGN with luminosity. This observational evidence implies 
that our simple models, where the torus does not evolve with 
neither luminosity nor redshift, might not be an accurate repre- 
sentation. 

Treister & Urry ({2005 ) found a simple explanation for 
Barger's results assuming a Modified Unification Model, where 
the percentage of obscured AGNs varies linearly from 100% at 
L x = 10 42 erg sec" 1 to 0% at L x = 3 ■ 10 46 erg sec" 1 . We can ex- 



amine how this luminosity evolution would change our results 
for the torus models. Using the luminosity dependency for the 
fraction of Type II AGN determined by La Franca et al. ((2005), 
we find that in the case of Model 1 (the one with the best sta- 
tistical results), parameter y, which is mostly responsible for 
the variations observed in this fraction, changes to m 7.0, 14 
and 23 (fixing R m /R t = 1.1) for the luminosity bins 10 42 " 43 , 
10 43-44 and 10 44 " 45 erg sec" 1 , respectively. These values take 
into account the observational incompleteness for the Type II 
population in La Franca's X-Ray sample (~l/3 of the sources 
would not observed). Table 1 shows that our results give an 
estimate of 7.0 < y < 9.0, roughly coinciding with y for the 
lowest and middle luminosity bins in La Franca's sample. 

Unfortunately, the current statistics of the observed Nh dis- 
tribution does not allow us to test such luminosity dependency. 
However, the number density of the most luminous sources is 
much smaller than those of low and intermediate luminosity, 
with fractions of 63%, 30%, and 7% for the 10 42 " 43 , 10 43 " 44 and 
10 44 " 45 erg sec" 1 luminosity bins (integrating Lx from 10 42 to 
10 45 erg sec" 1 in the Ueda's luminosity function). Therefore, 
the general Nh distribution is clearly dominated for the less 
powerful (Seyfert like) AGN, for which the derived torus prop- 
erties might be a good representation. 

5.2. Column Densities Uncertainties 

Our analysis is based on the distribution of the determined col- 
umn densities. From this measurements, it is not possible to 
determine how far from the central source the absorption is 
produced. This implies that overdensities of gas located closer 
than the graphite sublimation radius will not have the assumed 
Gas-to-Dust ratio. This possibility is consistent with the ex- 
treme variability of the column density observed in the Seyfert 
galaxy NGC 1365 (Risaliti et al. 2005). Also, assuming higher 
metallicities gives a harder photoelectric absorption, increasing 
the estimation of the number of Type I sources. Finally, chang- 
ing the parameter R y from 3 . 1 to 5 .0 produces small differences 
in our results. 

A spectroscopic analysis such as the one conducted by 
Tozzi et al. (120061 1 is essential to estimate the distribution of 
column densities since particular spectral features for each 
source can be present in the data. Broad band estimations, us- 
ing Hardness Ratios are not precise enough and are highly 
affected by the template assumed for the X-Ray source. These 
uncertainties affect the best-fit model parameters, and the er- 
ror bars for these estimations. We conclude that a direct fit to 
the spectral features is determinant at the time to describe the 
obscuring properties of the AGN population. 

5.3. Torus Model Results 

Recent IR observations have derived compact torus structures 
in nearby Seyfert II galaxies (Jaffe et al. 2004; Prieto et al. 
2004 2005 1 Swain et al. 2003 ) in agreement with the sizes pre- 
dicted for model 1, while for models 3 and 4 this physical scale 
can be freely adjusted. 
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It has been argued that a constant torus density distribution 
does not explain the reemission IR spectrum in AGN, but the 
inclusion of a clumpy medium solves this problem. Models 1 
to 3 can be discretized using Equation [1] while adopting small 
optical depths per cloud allows predictions of Type I sources at 
^ 10 21 cm" 2 , avoiding a large dichotomy between Type II and 
Type I objects (see Model 4 for details). 

Exponential angular dependencies of the parameter p (see 
Equation [TJ allow smooth density decrements from the equa- 
torial plane to the poles, so that a wide range of Nh can be 
obtained, implying that Model 1 and 3 are the most favorable 
obscuring structures. On the other hand, radial dependencies 
for the models (J3 + 0) can shrink the effective torus radius, in- 
creasing the inner equatorial density, more according with dy- 
namical simulations of a thick clumpy torus (see Beckert & 
Duschl|2004]). 

5.4. Highly absorbed sources 

According to the X-Ray analysis of Seyfert II galaxies pre- 
sented by Maiolino et al. (119981 1, the number of Compton Thick 
sources represents a large fraction of the AGN population. In 
this work, the assumption of a maximum column density at 
Nh = 10 25 cm" 2 allows us to estimate this highly absorbed 
population from our geometrical assumptions (see Figure |4j. 
We roughly predict intrinsic Compton Thick source fractions 
of 27%, 3%, 26% and 58% for Model 1, 2, 3 and 4 respec- 
tively. Work by Tozzi et al. (2006) estimate a lower limit for 
the Compton Thick fraction of ~ 10% in the CDF-S, while 
Beckmann et al. (2006) find 4 Compton Thick sources of a to- 
tal of 40 AGN observed with INTEGRAL in the 20-40 keV 
energy range at z ~ 0. These values are somewhat smaller, but 
still consistent with predictions from Models 1 and 3 given the 
uncertainties involved. 

6. Conclusions 

1 . Using the observed distribution of Nh columns for AGNs 
in the CDF-S, we have constrained the geometry of the ob- 
scuring region present in AGNs. Four different torus ge- 
ometries, based on previously proposed models (see Figure 
[T), have been explored. No luminosity or redshift depen- 
dency has been included in the torus modelling. 

2. Results from our fit to the Nh distribution are presented in 
Table Q] which shows the model parameter best-fit values. 
These results can be used to estimate physical properties of 
the torus, such as the inner and outer radius, and the density 
at the equatorial plane, as a function of the AGN luminosity. 

3. Detailed analysis shows that a matter density profile (p(4>)) 
with an exponential dependency in 0, such as in Model 1 
and 3, gives the best representation of the wide range of 
observed column densities (from ~ 10 20 cm" 2 to ~ 10 25 
cm -2 ). A constant density distribution (Model 2) requires 
very extended structures to recover the wide observed Nh 
range (since Nh 00 As), not according with actual obser- 
vations of local Seyfert galaxies where pc-scale torus-like 
structures have been detected. It also underestimates the 
fraction of Type II sources. A Gaussian angular dependency 



(Model 4), as in the clumpy model proposed by Nenkova et 
al. (|2002), results in a substantial overestimation of highly 
obscured sources. We also find that the optical depth for in- 
dividual clumps, Tvu has to be ~ 60 times smaller than the 
value adopted by Nenkova to compute the torus IR reemis- 
sion. We note that the clumpy torus treatment could also be 
applied to other model distributions using Equation Q] and a 
suitable value for -ry,-. Models 3 and 4 suffer of multiple x 1 
minima caused by the dichotomy between Type I and Type 
II sources (given by the critical value of <p and the presence 
of at least 1 clump along the line of sight, respectively). 
4. We therefore conclude that Model 1, a classical 'donut 
shape' , is the best parameterization for the obscuring region 
around AGN. Since no luminosity dependency has been in- 
cluded in our modelling, the results can be regarded as a 
luminosity averaged parameterization of the torus, which 
is largely dominated by low and intermediate luminosity 
sources. 
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Fig. 1. Geometrical matter density distributions assumed for the torus models. Figures a), b) c) and d) are based on the previous 
work by Treister et al. J2004I I. Pier & Krolik ( U992K Granato & Danese ( 119941 1 and Nenkova et al. (2002 ), respectively. 




Fig. 2. Best-fit parameters found for each model following Section [4] Contours represent 1, 2 and 3 <x confidence limits. Note 
that for Model 1 = 0.0. 
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Torus Model 1 
Intrinsic NH distribution' (Tozzi et al. 2006)' 
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Fig. 3. Comparison between the column densities of the observed Chandra sources in the CDF-S field (crosses) and the best-fit 
results obtained from our modelling. For each model, two histograms are presented varying the range of interest of one parameter 
only. 
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Fig. 4. Geometrical matter density distribution using best fit parameters for each model. Model 1: R m /R t = 1.14, y = 9.0 and 
/3 = 0.0; Model 2: hi a = 3.0 and b/a = 350; Model 3: <f> c = 1.1 and y = 9.5; Model 4: cr = 0.4 and r vi « 1 (continuous case). 
The inner radius is set at 1 .2 pc for all models, and an outer radius at 12 pc for models 3 and 4 (but adjustable for more compact 
structures). Different shades mean Alog(A^ H ) = 0.5, starting from 10 25 crrT 2 at the equatorial plane. The dashed line represents 
the line of sight that divides the Type I and the Type II sources. Note that Model 2 is shown at a different scale. 



